Dissipation Mechanisms in Thermomechanically Driven Silicon Nitride Nanostrings 
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High-stress silicon nitride nanostrings are a promising system for sensing applications because of their ultra- 
high mechanical quality factors (Qs). By performing thermomechanical calibration across multiple vibrational 
modes, we are able to assess the roles of the various dissipation mechanisms in these devices. Specifically, 
we possess a set of nanostrings in which all measured modes fall upon a single curve of peak displacement 
versus frequency. This allows us to rule out bulk bending and intrinsic loss mechanisms as dominant sources of 
dissipation and to conclude that the most significant contribution to dissipation in high-stress nanostrings occurs 
at the anchor points. 



The extremely high values of mechanical Q that have been 
reported in silicon nitride nanostrings 1 ™ have generated a 
great deal of excitement in the nanomechanics community. 4 
These devices are ideal for use in mass sensing, 5 temperature 
sensing, 6 and optomechanicP" 9 -! and have proved to be an in- 
valuable platform for research into the quantum properties of 
nanoscale resonators. 10 11 Nanostrings possess all the desired 
properties for these endeavors, including small mass, high fre- 
quency, and high Q, with correspondingly large displacement 
amplitudes. This combination makes them sensitive to exter- 
nal perturbation yet still within the limits of current detection 
techniques.^! 

In this manuscript, we demonstrate thermomechanically 
limited detection of up to six mechanical nanostring modes. 
Furthermore, we present a set of devices in which all harmon- 
ics fall upon a single curve of calibrated peak displacement 
versus frequency. As we discuss below, we are able to infer 
that the mechanical Q is limited by dissipation processes op- 
erating in the vicinity of the anchor points — thus suggesting 
ways to further engineer the mechanical Q. 

Silicon nitride nanostrings are devices under extremely 
high tensile stress (a = 0.8 GPa for our devices). The ac- 
cepted understanding is that the tension along their length re- 
sults in large stored elastic energy and a correspondingly large 
mechanical Q? Experiments have confirmed that Q increases 
with mechanical tensioning of non-prestr essed (low tension, 
non-stoichiometric silicon nitride) devices] 13 ! 14 ! corroborating 
the tension's central role. 

In addition to causing the high Q, the tension com- 
pels the devices to behave like strings rather than doubly 
clamped beams, as one would usually expect for this geom- 
etry. The harmonics are at integer multiples of the funda- 
mental frequency, 15 v n = nv\, where v\ is the fundamen- 
tal mode frequency and n indicates the mode number; the 
mode frequency depends only on the length of the string (not 
on the width or thickness). These features are in contrast to 
the more complicated harmonics of doubly clamped beams, 16 
which are realized in low-stress silicon nitride devices of the 
same geometry.^ Furthermore, it has now been shown that 
very high mechanical Qs exist for nanostrings under tension 
in a variety of other materials, including the polymer SU-8JSI 
aluminum, 6 19 and AuPd. 20 This demonstrates that the high Q 
is a byproduct of the tension and string-like behavior and not a 
special material property of the silicon nitride itself. It is likely 



that most materials that can withstand high tensile stress can 
be coerced toward higher Q in this way. 

Nanostrings offer significant potential for sensing 
applications, 21 provided that they can be accurately 
calibrated!^! In this work, the absolute displacement of 
a nanostring is determined from measurements of its thermal 
activationP 3 -! This procedure is often restricted to the funda- 
mental mode, since it is the mode of largest displacement. 
Because of the high Q of nanostrings, however, we are able to 
detect the thermomechanical motion of up to six harmonics 
for our longest strings. 

A key insight is that such multimode thermomechanical 
calibration reveals important information about whether the 
dominant dissipation processes in nanostrings are localized at 
the anchor points or in the bulk. The calibrated peak displace- 
ments of our devices fit an analytical form that follows from 
very general considerations of the leading order damping al- 
lowed by symmetry. Nanoscale devices may dissipate quite 
differently from macroscopic oscillators, 24 since these devices 
have periods of oscillation that are comparable to their thermal 
transport times and phonon mean free paths that rival their lin- 
ear size. But our analysis is agnostic with respect to the under- 




FIG. 1. (a) Schematic of our measurement set-up. The silicon 
nitride nanostring is mounted in an optical-access vacuum chamber, 
(b-g) Cartoons showing the mode shapes of the first six harmonics 
of a nanostring. (h-i) Scanning electron micrographs of the 215 /im 
long silicon nitride nanostring, with scale bars of length 20 fim and 
3 fim respectively. 
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FIG. 2. Calibrated thermomechanical spectra of the first six har- 
monics of the 215 fim long nanostring (the device shown in Fig.[TJ. 
These are single sweeps and are post-processed by a sliding average 
over five points. The data are fit to an analytical expression for the 
power spectral density function (see Ref. 29.) 



lying microscopic mechanism and is consistent with a number 
of different scenarios for losses at t he a nchor points (includ- 
ing dissipation by phonon tunneling 25 26 and other models of 
local bending at the anchor pointsP^) 

Our silicon nitride nanostrings are fabricated from stoichio- 
metric silicon nitride deposited onto silicon dioxide on a sili- 
con handle using LPCVD. This is cost effective as wafers are 
commercially available for « $50 for a 100 mm wafer (Rogue 
Valley Microdevices). We perform standard optical lithog- 
raphy followed by reactive ion etching. The nanostrings are 
released by a buffered oxide etch of the 2 /im thick silicon 
dioxide layer. 

The completed devices, shown in Figs. 1(h) and l(i), are 
mounted onto a piezoelectric, which can provide mechani- 
cal actuation when excited by the lock-in amplifier (Zurich 
Instruments HF2LI). The experiment is housed in a short- 
working-distance vacuum chamber at s» 2 x 10~ 6 torr, which 
is sufficiently low to eliminate pressure related damping. 14 
The displacement of the nanostring is measured via opti- 
cal interference with respect to the silicon substrate (HeNe 
632.8 nm focused to an ~ 1 pm diameter spot), as shown 
in Fig. l(a)P3The 100 pW of optical power impinging on 
the string causes no significant heating: we have verified that 
the resonance frequencies do not shift as a function of opti- 
cal powerf-^ a consequence of the low absorption of silicon 
nitride in the visible. Interferometric detection can be per- 
formed with or without a voltage applied to the piezo. In the 
latter case — appropriate for measurement of the thermome- 
chanical actuation — the HF2LI acts as a spectrum analyzer 
with a user-controlled measurement bandwidth. Because the 
Qs of the devices are so high, a bandwidth of 1 Hz is required. 
We have identified just how narrow the bandwidth must be by 
ensuring that the Q measured via thermomechanical actuation 
is consistent with both the Q of the driven device and the Q 
extracted from a ringdown measurement versus time. 

In a simple model, the string is described by its vertical 



displacement z(x, t), which is a function of the horizontal po- 
sition x along the string and the time t. The string endpoints 
are fixed at z(0, t) = z(L, t) = 0, and the string achieves its 
natural length L when the displacement of the string is every- 
where zero. In the linear regime — applicable for all measure- 
ments reported here — the length of the string deviates from 
L by an amount f^dxz 2 .. (Here z x — dz(x,t)/dx denotes 
the partial derivative.) By virtue of the strong intrinsic tensile 
stress a of the silicon nitride, there is a corresponding restor- 
ing force proportional to the local curvature when the string is 
stretched. Hence, the string obeys a modified wave equation 
pz u = crz xx — D[z(x, t)}, where the functional D encodes the 
various dissipative processes in operation. 

It is convenient to decompose the motion of the system into 
its normal modes 



z(x, t) = a n (t) sin 



/rnrx\ 



(D 



which are indexed by a positive integer n, the mode number. 
The first six mode shapes are sketched in Fig. l(b)-(g). Each 
mode acts as a nearly independent damped harmonic oscilla- 
tor obeying 



d 2 < 



dt 2 



+ ln-^j- + muj n a n = /„(*). 



(2) 




Here, m is the geometric mass of the string, u) n the angular 
frequency, and j n the damping coefficient (7„ = mu n /Q n ). 
In the case of thermomechanical actuation, the forcing term 
f n (t) can be understood as a continuous-time stochastic pro- 
cess that acts as a source of thermal noise. The sequence of 
resonance frequencies 



(3) 



depends on the device length and on the intrinsic tensile stress 
and density p of the silicon nitride. 8 Values in the literature for 
the density vary, and therefore we use our measured frequen- 
cies and lengths, along with the tensile stress reported by the 
manufacturer, a = 0.8 GPa, to determine the density for our- 
selves. Our estimate of 3000 kg/m 3 is consistent with other 
reported values 

We now briefly discuss the form of the damping coefficient. 
Since the experimental chamber is evacuated to m 2xl0 -6 torr, 
the dissipation cannot have a substantial contribution from 
any viscous term z t that depends on the absolute motion of 
the string.^ Instead, the leading order contribution from in- 
ternal material processes must go as z xx ,u which describes 
how fast a volume element in the string is moving with re- 
spect to neighboring ones. The exception is in the vicinity of 
the anchor points, where dissipation can depend on the time 
rate of change of the angle that the string makes with respect 
to its connection point (z x ,t)- By this reasoning, the effective 
Q for each mode is well approximated by 
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Here, we have taken 7„ to be a sum of three contributions 
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damping, damping at the anchor points, and damping due to 
bending of the bulk string. (The tilde-decorated quantities are 
material constants that have had their leading order u n depen- 
dence factored out.) Note that each contribution appears in the 
denominator of Eq. Q with a characteristic mode number (or, 
equivalently, resonance frequency) dependence, whose expo- 
nent p — 1 tracks the number of x and t derivatives in the 
relevant dissipation term d p+1 z/dx p dt. In particular, Q val- 
ues that are roughly constant across modes are indicative of 
dissipation at the anchor points; dissipation originating deep 
in the string, on the other hand, should lead to Q that decays 
inversely with the resonant frequency. The former behavior 
is consistent with our understanding that the relative contribu- 
tion ^buiky- anchor ^ (j-i/2 j s lively to be small in high stress 
materials. 

Thermomechanical calibration allows us to quantitatively 
compare these loss channels. Once a thermally driven spec- 
trum is acquired, as in Fig. 2, we calibrate the voltage signal 
to an absolute displacement via the equipartition theorem. Fit- 
ting the power spectral density function 29 to the square of the 
voltage signal divided by the measurement bandwidth, we ex- 
tract four parameters: the displacement noise floor, the Q, the 
resonance frequency, and the factor a (having units m 2 /V 2 ) 
that converts between the voltage signal and the displacement. 
The only uncertainty is the mass of the nanostring, which we 
compute as m = pLwt. The thickness t = 250 nm is the 
same in all devices, whereas the width w is idiosyncratic. Ac- 
curate knowledge of the density and geometry is important for 
this calibration. The main results are the peak displacement of 
the beam at the resonance frequency (Fig. 3a), and the qual- 
ity factor Q n (Fig. 3b), both of which are extracted from the 
fit to the calibrated spectrum. Calibration is repeated for each 
measurement of a resonance spectrum, since a varies with the 
location of our laser spot. 

The peak displacements and frequencies of the second 
mode of the 215 /im device and the first mode of the 102 /im 
device are nearly identical [see Fig. 3(a)]. The coincidence 
of the frequencies (following from the near identity 1/102 pa 
2/215) is a consequence of the n/L scaling in Eq. ([3]). As is 
characteristic of strings, the frequency is proportional to the 
sound velocity and mode number but otherwise depends only 
on the string length and not on any of its transverse dimen- 
sions. (The width of the nanostring has been shown to be 
important when w > 3 /tmHlbut when 55 nm < w < 1.5 /an 
the Q is independent of w at low pressures. 14 ) There are addi- 
tional degeneracies — the 3 rd mode of the 215 /mi device and 
the 1 st mode of the 74 /im device; the 4 th mode of the 215 /tm 
device and the 2 nd of the 102 /im device; the 6 th mode of the 
215 /im device, the 3 rd of the 102 /tm device, and the 2 nd of 
the 74 /tm device — that together constitute proof of these sys- 
tems' string-like behavior. 

On the other hand, it is not obvious that modes with identi- 
cal frequencies should also have identical peak displacements, 
as seen in Fig. 3(a). We note this behavior of the peak dis- 
placement is not universal. Rather, it is simply fortuitous that 
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FIG. 3. (a) Peak displacement and (b) Q from thermomechanical 
calibration of three devices, all 250 nm thick. Black circles are for a 
string 215 /tm long and 2. 1 /im wide, green squares for 102 /im long 
and 1.1 /im wide, red triangles for 74 /im long and 1.2 /im wide. The 
blue curve in (a) is a fit to Eq. 6. 



we have devices that fall on the same curve, and we have 
measured other devices that do not fit this trend (for exam- 
ple devices measured at higher pressures). The fact that these 
devices were manufactured identically, from the same wafer, 
and measured under rigorously controlled conditions, aids in 
the coincidence. In Fig. 3(b), the harmonics of a longer device 
are seen to have higher Q (lower dissipation) than the modes 
of a shorter device with identical frequencies. Conveniently, 
there exist in our data set modes from different devices that 
have nearly identical frequencies and peak displacements but 
substantially different Qs. We must conclude that two such 
modes differ only in losses at the anchor points. Specifically, 
since they have identical peak displacements — which is only 
the case because they are thermomechanically driven — these 
modes have identical local curvature, and therefore the Q can- 
not be limited by dissipation due to bulk bendingP^l We would 
be unable to make such a claim for a general actuation sce- 
nario in which the Q does not depend on amplitude, since both 
the stored energy and dissipation increase with driven ampli- 
tude. In our thermomechanical scenario, however, the only 
difference between such modes is the ratio of anchor points 
(two) to the number of lossless nodes (n — 1) [see inset to 
Fig. 3(a)]. The smaller this ratio, the smaller the effective 
loss at the anchors — dissipation is dominated by the anchor 
points, be it through local curvature^^ or the tunneling of 
phononsPSD 

Moreover, our analytical results tell us that the peak dis- 
placement versus frequency curve contains quantitative infor- 
mation about the dissipation! 2 ^ In particular, if it is true that 
loss at the anchor points dominates, then the peak displace- 
ment should decay as A more precise statement of 
the predicted frequency dependence is as follows: the spectral 
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density of the squared displacement is 

the corresponding displacement Az = ^ S n {i>™ x )Ai> over 
an arbitrary frequency interval Av goes as 

/q^ 1 

V ^/Vt^'V 1 + (7 buUc /7 anchor )27r^ (6) 

1 

~ 3/2 . 5/2' 

The fit shown in Fig. 3(a) is of this form, with a value of r = 

TT^bulk^anchor _ (_q 3 ± 3.5) x 1Q-7 s tha{ j s effectively 

zero within measurement error. This limit, r _1 ^> v n for 
all measured n, verifies that the ratio of bulk bending losses 
to losses at the anchor points is negligible. This is consistent 
with the weak variation of Q with respect to mode number that 
we see in Fig. 3b. (When displacements are large, as they are 
for the 215 nm string, there is an additional even/odd effect 
related to the parity of the mode shape; this has been shown in 
high-stress silicon nitride membranes^ — the 2D equivalent of 
strings — to result from phonon tunneling through the anchor 
points.) Finally, we note that the rate of energy loss is strongly 
influenced by the design of the junction between the string 
and the support. 3 31 Fitting the calibrated peak displacement 
versus frequency of the thermomechanical motion, as we have 



done here, can now be considered a straightforward test of 
the dominant dissipation and can be applied to other anchor 
designs. 

Establishing that nanostrings under tensile stress over- 
whelmingly dissipate from their anchor points opens a path 
for further increasing Q. Specifically, one could design anchor 
points at the ends of the beams that either carefully suspend 
the device^ or reflect phonons back into the string instead 
of allowing them to leak into the substrate. This has been 
achieved in part in Ref. [3] by clever fabrication of the nanos- 
trings and sample mounting. A scheme to eliminate phonon 
tunneling loss at the anchor points would be to design anchor 
points with complete phononic bandgaps at the frequencies of 
the mechanical modes of the nanostrings. This was rec ently 
accomplished for very high frequency modes (> 3 GHz)P323 

In conclusion, we have fabricated high Q silicon nitride 
nanostrings that enable multimode thermomechanical calibra- 
tion. Measurements of the peak displacement versus fre- 
quency for a collection of devices collapse onto a single curve, 
even though the mechanical Qs of these devices do not. An 
analysis based on the leading order dissipation allowed by 
symmetry rules out the possibility that bulk bending and in- 
trinsic loss mechanisms play a leading role. Instead, we 
confirm that loss mechanisms at the anchor points dominate. 
Hence, careful engineering of the anchor points in high stress 
nanostrings may offer a path to more sensitive devices and 
new applications. 
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Appendix A: Supplementary Information 

The string under consideration has a uniform (volumetric) 
mass density p and a uniform cross-sectional area A. Even 
at rest and unstretched, the string feels a tension a A along its 
length because of a stress a that is intrinsic to the material. 

Any smooth string profile z(x) between x = and x = L 
traces out a total length 



\fda 



dz 2 



= f dxy/l 
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(d xZ y - -{d xZ f + 
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The string acquires excess potential energy whenever it is 
stretched beyond its natural length L. For small amplitude os- 
cillations, and in the absence of any explicit dissipation mech- 
anism, the energy stored in the string is 

r-L 



H 



dx 
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2p~A' 



-n 2 



Aa 



(d x z) 



(A2) 



Here H(x) is the momentum conjugate to the displacement 
field. The relevant boundary conditions are z(0) = z(L) = 0. 
The endpoints are not effectively clamped, so the gradients 
z x (0) and z x (L) are unrestricted. 

(To next leading order Eq. ( |A2| i contains a term propor- 
tional to EA(d x z) 2 , where E is the Young's modulus. In 
a high tensile stress system, the crossover from clamped- 
beam-like to string-like behave occurs over a small region 
AL ~ y/E/a < L. Thus, even when the fourth order term is 
non-negligible, the string picture is correct in a coarse-grained 
sense.) 

Variation of the Hamiltonian by way of the identity 

-(d x z(x)) 2 = -26(x - 
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leads to equations of motion 
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In other words, the string obeys <rd 2 z — pz. Written in terms 
of the normal modes, 
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the continuum wave equation reduces to a discrete (but infi- 
nite: n = 1, 2, . . .) set of ordinary differential equations 



p L 2 
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Each normal mode is wholly independent and executes har- 
monic motion with angular frequency w„ = (<r / ' pY^mr / L. 

We proceed under the assumption that there are processes 
that allow energy to be shared between all the modes (ergod- 
icity). Then, if the string is held in thermal equilibrium at 
temperature T = I/UbP, its properties are described by a 
partition function 



Z 



V[z(x)]V[IL(x)] e 
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which takes the form of a Boltzmann-weighted integration 
over the complete phase space of string shapes and momenta. 
The II(a;) contribution will not enter so long as we only mea- 
sure functions of the string displacement. For example, 
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Once again, we want to decompose the field into its nor- 
mal modes. This requires a formal change of variables in the 
partition function. We imagine imposing a resolution limit by 
breaking up the spatial domain into bins of width Ax = L/N. 
At each position x* = (j — |) Ax, the displacement field takes 

a value coarse-grained z(x*) = z 
The infinitesimals 
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are connected by a matrix J that itself has no dependence on 
{zj} or {a n }. Hence, the Jacobian of the transformation 
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yields only an inert constant that factors out of any expectation 
values. 

Substituting the normal mode expression into the potential 
energy term gives 
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El 2 



m 2 2 



(All) 
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Here, we've taken advantage of the orthogonality relation 

(A12) 



I 

Jo 



nirx rrnrx L 

dx cos cos = — d n . 

L L 2 ' 



Gaussian integration of the mode variables gives 



I[I\n' a n'] a n a m e 



(A13) 

which confirms that the thermally averaged potential energy 



Jo 



dx 



Ao 



(d x zf ) = 



Nk B T 



(A14) 



has k B T/2 per mode. 

There's an important subtlety here: the expectation value 
( a n) = seems to imply that ^uKa 2 ) = k B T is the 

potential energy in the n th oscillator. It is easier to make the 
connection with the conventional invocation of the equiparti- 
tion theorem (half k B T per degree of freedom) if we write 

1 m 1 1 

-^ 2 n (a 2 n ) = -ptJ 2 n (a 2 n ) = -k B T (A15) 

in terms of the properly normalized mode amplitude a n — 
(2/L) 1 / 2 a n (with awkward units of distance to the half 
power). Here, p = m/L is the linear mass density. 

The average squared displacement of the thermally agitated 
string is 



(z(xf) 



. nirx . rrnrx , . 
sin — — sin — — {a n a m ) 



n.m 

E 



2k B T 



rnrx 



(A16) 



■ sin 



mujt, 



Since the spectrum is harmonic (cj n = nui), it is straightfor- 
ward to evaluate this expression at the center of the string: 



(z(L/2) 2 ) = 



2k B T 



muz 



n 2 



(A17) 



_ 2k B T pL 2 ir 2 
~ pL ^ x Y 
_ k B TL 

~ 4(7 ' 

We are interested, however, is resolving the motion in 
the frequency domain. To do this, we formally identify the 
thermal and time averaged values (a„)therm = ( a ^)time an d 
make use of Parseval's theorem to convert to Fourier space, 
treating the previous result as sum rule on the total square- 
displacement spectral density: 



1 



to 



PC 

K(t)\ 2 = / 

Jo 



dui S n (uj) = 



2k B T 



mujt, 



(A18) 

For the undamped oscillator case we have considered so far, 

(A19) 



but we expect the peak to be broadened once dissipation is 
introduced. We can solve the relevant damped harmonic os- 
cillator 



~r . ■> _ f(t) 



On + 7T a « + UJ " a 
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(A20) 



[appearing as Eq. (2)] by applying the usual trick of driving it 
at a single frequency and complexifying the variable a n (t) — 
Rea n e luJt . From 



(-w + wja„ = — — = C 



we get 



C sin(a;»t + <j> n ) 

N /( w 2_ w 2 ) 2 + (wWri/Qn) 2 



(A21) 



(A22) 



with a phase shift <j> n — arctan[2a;a;„Q^ 1 /(a; 2 — uj 2 )]. The 
unknown C is constant (independent of ui) because the ther- 
mal noise is white. Its value can be fixed by taking the time- 
averaged square displacement 



(a„) 2 = 



\C 2 



(wl - UJ 2 ) 2 + (oJUJ n /Q n ) 2 

and equating the integral 



(A23) 



L 



dw 



\C 2 



= \c 2 ^k (A24) 



(w 2 -w 2 ) 2 + (ojuj n /Q n ) 2 2 2ujI 



with the known weight under the curve, 2k B T/muj 2 . This 
establishes that 



and in turn yields 

S n (u>) 



i g2 ^ 4fc B r 

2 mirQn 



Au n k B T 



mirQ n [(ujl - uj 2 ) 2 + (ujoj n /Q n ) 



21 ' 



(A25) 



(A26) 



The corresponding expression in terms of conventional fre- 
quency is 



S n (u) = 2vtS„H = 



v n k B T 



rmr*Q n [{u 2 n -v 2 ) 2 + {uv n /Q n ) 2 X 

(A27) 

At the peak of the distribution 



the value is 



s»(C) 




mvT 3 (l- 4^tW 



(A28) 



(A29) 



In the experimental situation, the relevant quantity is the 
power spectral density, obtained by squaring the measured 



7 



voltage near resonance by the bandwidth of the lock-in ampli- 
fier. It is related to the square-displacement spectral density 
by 



V 2 



Sy = = (noise floor) + aS z {v) 



«o + 



ak B T \ Vn 

mir* )QnM-v 2 ) 2 + (w n /Qn) 2 } 



(A30) 



where a is a conversion factor having units V 2 /m 2 . It is 
straightforward to extract a and s from a fit to the data. 
(There is a unique a associated with each peak.) Hence, from 
the measured voltage we report a root square-displacement 
spectral density according to 



V 



\fuKv 



(A31) 



